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Abstract 

We propose a method for computing the cohomology ring of three-dimensional 
(3D) digital binary-valued pictures. We obtain the cohomology ring of a 3D 
digital binary-valued picture J, via a simplicial complex K{I) topologically 
representing (up to isomorphisms of pictures) the picture /. The usefulness of 
a simplicial description of the "digital" cohomology ring of 3D digital binary- 
valued pictures is tested by means of a small program visualizing the different 
steps of the method. Some examples concerning topological thinning, the vi- 
sualization of representative (co)cycles of (co)homology generators and the 
computation of the cup product on the cohomology of simple pictures are 
showed. 

Keywords: Digital topology, chain complexes, cohomology ring. 
1. Introduction 

The homology groups (given in terms of number of connected components, 
holes and cavities in the 3D digital picture), the digital Euler characteristic or 
the digital fundamental group are well-known operations in Digital Topology 



16|,lll|. All of them can be considered as translations into the discrete setting 
of classical continuous topological invariants. In order to prove that a digital 
topology operation vr^, (associated with a continuous operation ttc) correctly 
reflects the topology of digital pictures considered as Euclidean spaces, the 
main idea is to associate a "continuous analog" C (/) with the digital picture 
I. In most cases, each digital picture / is associated with a polyhedron 
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C(J) [nl, [lil Eol, U. It is clear that C(J) "fills the gaps" between the black 
points of / in a way that strongly depends on the grid and the adjacency 
relations chosen for the digital picture /. Recent attempts to enrich the list 
of computable digital topological invariants in such a way can be found in 

0. 

In this paper, starting from a 3D digital binary-valued picture /, a sim- 
plicial complex K{I) associated with / is constructed, in such a way that 
an isomorphism of pictures is equivalent to a simplicial homeomorphism of 
the corresponding simplicial complexes. Therefore, we are able to define the 
digital cohomology ring of I with coefficients in a commutative ring G, as the 
classical cohomology ring of K{I) with coefficients in G (see [HI]). In order 
to compute this last algebraic object, it is crucial in our method to "connect" 
the chain complex C{K{I)) canonically associated to K[I) and its homol- 



ogy H{K{I)), via an special chain equivalence [15|: a chain contraction [14 
We will obtain this goal in several steps. Using the technique of simplicial 
collapses joj, we topologically thin K{I), obtaining a smaller simplicial com- 
plex Mt^pK{I) (with the same homology as K{I)) and a chain contraction 
connecting their respective chain complexes C{K{I)) and C{Mt^pK{I)). The 
following step is the construction of a chain contraction from C{MtapK{I)) 
to its homology H{MtapK{I)). Having all this information at hand, it is easy 
to compute the digital cohomology ring of I for a given commutative ring G. 
In this way, cohomology rings are computable topological invariants which 
can be used for "topologically" classifying (up to isomorphisms of pictures) 
and distinguishing (up to cohomology ring level) 3D digital binary-valued 
pictures. 

A small program for visualizing these cohomology aspects in the case 
G = Z/Z2 has been designed by the authors and developed by other^. 
This software allows us to test in some simple examples the potentiality and 
topological acuity of the method. 

We deal with digital pictures derived from a tessellation of three-space by 
truncated octahedra. This is equivalent to using a body-centered-cubic-grid 
whose grid points are the points (x, y, z) G in which x = y = c (mod 2) 
(see [l3|). The only Voronoi adjacency relation on this grid is 14-adjacency. 
Using this adjacency, it is straightforward to associate to a digital picture 



^ The 1st version was programmed by J.M. Berrio, F. Leal and M.M. Maraver [3|; the 



2nd version by F.Leal. ,http://www. us.es/gtocoma/editcup.zip 
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/, a unique simplicial complex K{I) (up to isomorphisms of pictures) with 
the same topological information as /. This grid is important in medical 
imaging applications due to its outstanding topological properties and its 
higher contents of symmetries. One advantage of the voxels in this grid is 
that they are more "sphere-like" than the cube, so that the volumetric data 
represented on this grid need fewer samples that on Cartesian cubic grid 17 . 

Since the objects considered in this paper are embedded in then the 
homology groups vanishes for dimensions greater than 3 and they are torsion- 
free for dimensions 0, 1 and 2 (see I, ch.lO]). The qth Betti number is defined 
as the rank of the gth homology group. In general, the 0th Betti number is 
the number of connected components, the 1st and 2nd Betti numbers have 
intuitive interpretations as the number of independent non-bounding loops 
and the number of independent non-bounding shells. According to the Uni- 
versal Coefficient Theorem for Homology, the Betti numbers are independent 



of the group of coefficients (see 15|, ch. 7]). Moreover, since the homology 



groups are torsion-free, the cohomology groups with coefficients in G are 
isomorphic to the homology groups with coefficients also in G (see [l^, ch. 
5]). Therefore, for simplicity we can consider that the ground ring is Z/Z2 
throughout the paper. Nevertheless, all the procedure we explain here, is 
valid for any commutative ring G. 

The paper is organized as follows. In Section 2, the technique associat- 
ing a simplicial complex to a 3D digital binary-valued picture is detailed. 
In Section 3, we explain a procedure for computing the cohomology ring 
of general simplicial complexes. In Section 4, we introduce the notion of 
digital cohomology ring of a 3D digital binary-valued picture and we show 
some examples concerning the visualization of representative (co) cycles of 
(co)homology generators and the computation of the cup product on the co- 
homology of simple pictures. Finally, Section 5 is devoted to conclusions and 
comments. 



2. Prom Digital Images to Simplicial Complexes 

Digital Images. We follow the terminology given in [l2| for representing 
digital pictures. A 3D digital binary-valued picture space (or, briefly, DPS) 
is a triple (V, /3, w), where V is the set of grid points in a 3D grid and each of 
/3 and w is a set of closed straight line segments joining pairs of points of V . 
The set /3 (resp. the set u) determines the neighbourhood relations between 
black points (resp. white points) in the grid. A 3D digital binary-valued 
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Figure 1: The 14-neighbours of a grid point p of the BCC grid. 

picture is a quadruple / = (V, /3, B), where (V, /3, u;) is a DPS and B (the 
set of black points) is a finite subset of V. 

An isomorphism of a DPS (Vi,/3i,a;i) to a DPS (V2,/32,i^2) is a homeo- 
morphism h of the Euclidean 3-space to itself such that h maps Vi onto V2, 
each /3i-adjacency onto a /32-adjacency and each wi-adjacency onto an U2- 
adjacency, and maps each /32-adjacency onto a /3i-adjacency and each 
U2 adjacency onto an wi-adjacency. An isomorphism of a picture Ii = 
(Vi, (3i,ijJi, Bi) to a picture I2 = (V2, /325 1^2, -B2) is an isomorphism of the 
DPS {Vi,(3i,uji) to the DPS {V2,/32,uj2) that maps Bi onto B2. 

The DPS used in this paper, is the 3D body-centered cubic grid (BCC 
grid) [l2|: The grid points V are the points (a, 6, c) G such that a = b = c 
(mod 2). The 14-neighbours of a grid point p with coordinates (a, 6, c) are: 
(a ± 2, 6, c), (a, 6 ± 2, c), (a, 6, c ± 2), (a ± 1, 6 ± 1, c ± 1). 
Simplicial Complexes. The four types of non-empty simplices in are: a 
0-simplex which is a vertex, a 1-simplex which is an edge, a 2-simplex which 
is a triangle and a 3-simplex which is a tetrahedron. In general, considering 
an ordering on a vertex set V, a g-simplex with g + 1 affinely independent 
vertices f < • ■ ■ < f g of is the convex hull of these points, denoted by 
{vq, . . . ,Vg). If i < g, an i-face of a is an i-simplex whose vertices are in the 
set {vq, . . . ,Vg}. A facet of a is a (g — l)-face of it. A simplex is shared if it is 
a face of more than one simplex. Otherwise, the simplex is free if it belongs 
to one higher dimensional simplex, and maximal if it does not belong to any. 



4 



A simplicial complex K is a collection of simplices such that every face 
of a simplex of K is in K and the intersection of any two simplices of is a 
face of each of them or empty. The set of all the g-simplices of K is denoted 
by K^'^\ A subset K' C K is a subcomplex of K if it is a simplicial complex 
itself. 

Let K and L be simplicial complexes and let \K\ and |L| be the subsets of 
R'* that are the union of simplices of K and L, respectively. Let / : K^^^ 

be a map such that whenever the vertices vq, . . . ,Vn oi K span a simplex 
of K, the points /(fo), • • • , f{vn) are vertices of a simplex of L. Then / can 
be extended to a continuous map g : \K\ — )■ \L\ such that if x = J2tiVi then 
g{x) = J2tif{vi). The map g is called a simplicial homeomorphism if / is 
bijective and the points /(t>o), . . . , /(fn) always span a simplex of L. 

Simplicial Representations. Given a 3D digital binary-valued picture 
/ = (V, 14, 14, B) on the BCC grid, there is a process to uniquely associate 
a 3-dimensional simplicial complex K{I). This simplicial complex is con- 
structed on the triangulation of the Euclidean 3-space determined by the 
previous 14-neighbourhood relation. The simplicial representation K{I) of 
the digital picture / is described as follows: consider the lexicographical or- 
dering on V (if f = (a, b, c) and w = (x, y, z) are two points of V, then 
V < w a a < X, oi a = X and b<y, or a = x,b = y and c < z). The 
vertices (or 0-simplices) of K{I) are the points of /. The i-simplices of K{I) 
{i G {1,2,3}) are constituted by the different sorted sets of i 14-neighbour 
black points of / (analogously, we could construct another simplicial complex 
whose i-simplices are the different sets of i 14-neighbour white points of /). 

Example 2.1. Consider the digital picture J = (V, 14, 14, i?) where B is 
the set {vo = (-l,-l,l),t;i = (-l,l,l),t;2 = (0,0,0),t;3 = (0,0,2),t;4 = 
(0, 2, 0)} ; then K{J) is the simplicial complex with set of maximal simplices 
{{vq, vi, V2, vs), {vi, V2, Vi)} (see FigureW\j. 

In the next section, we give a satisfactory algorithmic solution to the 
problem of the computation of the cohomology ring of finite simplicial com- 
plexes. This positive solution together with the naive simplicial construction 
described above will allow us to "cohomologically control" 3D digital binary- 
valued pictures (up to isomorphisms of pictures), since the following result 
holds. 

Theorem 2.2. Two digital binary-valued pictures, Ii = (V, 14, 14, i?i) and 
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Figure 2: On the left, the black points of the digital picture J and, on the right, the 
simplicial representation K{J). 



I2 = (V, 14, 14, B2), are isomorphic if and only if their simplicial represen- 
tations K{Ii) and K{l2) are simplicially homeomorphic. 

The proof of this theorem is straightforward and left to the reader. 

3. Computing the Cohomology Ring of Simplicial Complexes 

First of all, we briefly explain the main concepts from Algebraic Topology 



we use in this paper. Our terminology follows Munkres book [15[. In the 



next subsections, we reinterpret classical methods in Algebraic Topology in 



terms of chain contractions [IJ] that will enable us to design an algorithm 
for computing the cohomology ring of general simplicial complexes. 
Chains and Homology. Let i^' be a simplicial complex. A q-chain a is 
a formal sum of simplices of K^'^\ Since the group of coefficient is Z/Z2, 
a g-chain can be seen as a subset of g-simplices of K; the sum of two q- 
chains c and d is the symmetric difference of the two sets cU d and c (1 d. 
The g-chains form a group with respect to the component-wise addition 
mod 2; this group is the qth chain group of K, denoted by Cq{K). There 
is a chain group for every integer q > 0, but for a complex in R'^, only 
the ones for < g < 3 may be non-trivial. The boundary of a g-simplex 
a = {vq, ■ ■ ■ ,Vq) is the collection of all its facets which is a (g — l)-chain: 
= I]i=o('^0) ■ ■ ■ . . . ,Vq), where the hat means that Vi is omitted. 
By linearity, the boundary operator dq can be extended to g-chains. The 
collection of boundary operators connect the chain groups Cq{K) into the 
chain complex C{K): ■■■ % C^{K) % C2{K) % Ci{K) % Co{K) % 0. 
A g-chain a G Cq{K) is called a q-cycle if dq{a) = 0. If a = (9g+i(a') for 
some a' G Cqj^i{K) then a is called a q-boundary. We denote the groups of 
g-cycles and g-boundaries by Zq and Bq respectively. An essential property 
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of the boundary operators is that the boundary of every boundary is empty, 
dqdg^i = 0. This imphes that Bg C Zg for q > 0. Define the qth homology 
group to be the quotient group Zg/Bg, denoted by Hq{K). Given a e Zg, 
the coset a + Bg is the homology class of Hg{K) determined by a. We denote 
this class by [a]. For a complex K in R^, only Hq{K) for < g < 2 may be 
non-trivial. 

Cochains and Cohomology. With each simplicial complex K, we have 
associated a sequence of abelian groups called its homology groups. Now, we 
associate with K another sequence of abelian groups called its cohomology 
groups. They are geometrically much less natural than the homology groups. 
Their origins lie in algebra rather than in geometry; in a certain algebraic 
sense, they are "dual" to the homology groups. 

Let i^' be a simplicial complex. The group of q-cochains of K with 
coefficients in Z/Z2 is the group C^K) = {c : Cg{K) ^ Z/Z2 such 
that c is a homomorphism}. Observe that a g-cochain c can be defined 
on the g-simplices of K and it is naturally extended to Cg{K). There- 
fore, a g-cochain can be expressed as a formal sum of elementary cochains 
a* : Cg{K) — Z/Z2 whose value is 1 on the g-simplex a & K and on 
all other g-simplices of K. The boundary operator dg+i on Cg+i{K) in- 
duces the cohoundary operator 6g : C^i^K) C^^^i^K) via 5g{c) = cdg+i, 
so that Sg raises dimension by one. The collection of coboundary opera- 
tors connect the cochain groups C'^{K) into the cochain complex C*{K): 

C\K) % C\K) % C\K) h C'iK) % - --.We define to be the 

kernel of 6g and B'^^^{K) to be its image. These groups are called the group 
of q-cocycles and q-coboundaries, respectively. Noting that Sg = because 
92 = 0, define the gth cohomology group, H'^{K) = Z'i{K) / B'i{K) for g > 0. 

The cochain complex C*{K) is an algebra with the cup product C^{K) x 
C'i{K) CP+^iK) given by: 

(c ^ c){a) = c{{vo, ...,Vp))» c'{{vp, . . .,Vp+g)) 

where a = {vq, . . . , Vp+g) is a (p + g)-simplex and • is the natural product 



defined on Z/Z2 [15|, p. 292]. It induces an operation ^: Hp{K) x H'^iK) 
HP~^'^[K), via [c] ^ [c'] = [c c'], that is bilinear, associative, commutative 
(up to a sign if the ground ring is not Z/Z2), independent of the ordering of 
the vertices of K and topological invariant (more concretely, homotopy-type 
invariant) [lil, p. 289], since the coboundary formula 5p+g(c ^ c') = Sp{c) ^ 
c' + c^ Sg{c') holds for any c G Cp{K) and c' G C"?(ir). 
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Figure 3: A triangulation of a torus. 



Figure 4: On the left, the 1-cycles a and b, and on the right, the 1-cocycles c and d. 

Example 3.1. Consider the complex K pictured in Figure\^ which is ob- 
tained from a triangulation of a torus. 

It is easy to check that the two 1-chains a = (3,7) + (7,9) + (3,9) and 
b = (3,7) + (6,7) + (6,8) + (8,9) + (3,9) (see F%gure\^ are 1-cycles. For 
example, di{a) = (3) + (7) + (7) + (9) + (3) + (9) =0. Moreover, a and 
b are homologous: 92((6,7, 8) + (7,8,9)) = a + b. On the other hand, c = 
(2, 3)* + (3, 6)* + (6, 7)* + (7, 8)* + (8, 9)* + (2, 9)* and d = (4, 5)* + (5, 6)* + 
(6, 8)* + (7, 8)* + (7, 9)* + (4, 9)* are two 1-cocycles. To check this, we have to 
verify that Si{c) and 6i{d) vanishes on all the 2-simplices of K. For example 
(5i(c)((2,3,6)) = c(92((2,3,6))) = c((2,3)) + c((2,6)) + c((3,6)) = 0. To 
check that both 1-cocycles are not coboundaries is a more difficult task since 
we have to verify that S2{f) ^ c,d for any f being a 2-cochain. 

The cup product of c and d is a new 2-cocycle c-^ d. By direct computa- 
tion, we have that c d = (6,7,8)* . We obtain this by applying it on all the 
2-simplices of K. For example, (c ^ d){{6,7,8)) = c((6, 7)) •d{{7,8)) = 1 
and (c ^ d){{7, 8, 9)) = c((7, 8)) • d{{8, 9)) = 0. 

The example illustrates that while we can think of a 1-cycle as being a 
closed curve, the best way to think of a 1-cocycle is a picket fence. 
Chain Contractions. In a more general framework, a chain complex C is 

a sequence ■ ■ ■ C3 C2 Ci Co of abelian groups Cq 
and homomorphisms dq, indexed with the non-negative integers, such that 
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for all g, dgdg^i = . The qth homology group is the quotient group Ker 
dg/ Im dg+i, denoted by Hg{C). Let C = {Cg,dg} and C = {C'g,d'g} be two 
chain complexes. A chain map / : C — C is a family of homomorphisms 
{fq '■ Cq — )■ Cg}g>o such that d'gfq = fq-idg . A chain map f : C ^ C induces 
a homomorphism /* : H{C) — ?■ H{C'). 

Let us emphasize that a fundamental notion here is that of chain con- 
traction. 

Definition 3.1. Jl^ ] A chain contraction of a chain complex C to another 
chain complex C is a set of three homomorphisms (/, g, 0) such that: 

• f : C ^ C and g : C ^ C are chain maps. 

• fg is the identity map of C . 

• (f) : C ^ C is a chain homotopy of the identity map id^ of C to gf , that 
is, (pd + d(f) = idc + gf ■ 

Important properties of chain contractions are: 

• C has fewer or the same number of generators than C. 

• C and C have isomorphic homology groups [15, p. 73]. 

Let us recall that the key of our method for computing the cohomology 
ring of chain complexes, is the construction of chain contractions (/,(?, 0) 
of a given chain complex C to another chain complex "H (isomorphic to the 
homology of C). In this case, for each cycle a G C, the chain /(a) G H 
determines the homology class of a. Conversely, for each a G (which cor- 
responds to a homology class of H{C)), g{a) G C determines a representative 
cycle of it. Finally, if a G C is a boundary, then a' = (f)[a) is a chain in C 
such that d{a') = a. 

3.1. Topological Thinning 

Topological thinning is an important preprocessing operation in Image 
Processing. The aim is to shrink a digital picture to a smaller, simpler picture 
which retains a lot of the significant information of the original. Then, further 
processing or analysis can be performed on the shrunken picture. 

There is a well-known process for thinning a simplicial complex using 
simplicial collapses Q. Suppose K is a. simplicial complex, cr G i^T is a 
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maximal simplex and a' is a free facet of a. Then, K simplicially collapses 
onto K — {a', a}. An important property of this process is that there exists 
an explicit chain contraction of C{K) to C{K — {a', a}) joj. More generally, 
a simplicial collapse is any sequence of such operations. A thinned simplicial 
complex MtopK is a subcomplex of K with the condition that all the faces of 
the maximal simplices of Mt^^K are shared. Then, it is obvious that it is no 
longer possible to collapse. 

The following algorithm computes M^^^K (first step) and a chain contrac- 
tion (/top, Qtop, (ptop) of C{K) to C{Mt„pK) (second step). In particular, recall 
that this means that the (co)homology of K and MtopK are isomorphic. Each 
step of the algorithm runs in time at most 0{m?) if K has m simplices. 

Algorithm 3.2. Topological Thinning Algorithm. 
First step: Simplicial collapses. 

Input: A simplicial complex K. 
Initially, M^^pK := K, collapse := ( ), pair := True. 
While pair is True do 
pair := False. 
For each a G M^^pK do 

If a is maximal with a free facet a' in Mt^^K then 
M,^pK:=K-{a',a}, 
collapse := (cr', a) U collapse, 
pair := True. 
End if. 
End for. 
End while. 

Output: the simplicial complex: M^^^K 

and the sorted set of simplices: collapse. 

Second step: the computation of the chain contraction. 

Input: The simplicial complexes K and M^^^K 

and the sorted set collapse = (cr^, cxi, cr^, (j„). 
Initially, ft,p{a) := a, 0t„p((T) := for each a e K ; 

and Qtopio') := a for each a G Mt^pK . 
For i = 1 to i = n do 

ftop{(T'i) := ftop{dai + a[), 
0top(o-D := CTi + (ptcpidcTi + 0-^), 
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Figure 5: The simplicial complexes L (on the left) and M^^pL (on the right). 

ftopicTi) ■■= 0. 
End for. 

Output: the chain contraction {Up, 9topAtop) of C{K) to C{M,,^K). 

Example 3.3. Consider the simplicial complex L whose set of maximal sim- 
plices is {(1, 5), (2, 5), (1, 2, 3), (2, 3, 4)} (see Figure\^. Applying the first 
part of the algorithm above we have that Mt^^L = {(1,3), (3,4), (2,4), (1,5), 
(2,5)} anc? collapse = ((2,3), (2,3,4), (1,2), (1,2,3)). The stages of the sec- 
ond part of the algorithm is showed in the following table: 





ftopitr) 


^top{cr) 


(2,3) 


/*..(2,4) + /,„,(3,4) 
= (2,4) + (3,4) 


(2,3,4) + (^,„,(2,4) + (/.,„,(3,4) 
= (2,3,4) 


(2,3,4) 








(1,2) 


/*.p(l,3) + /..,(2,3) 
= (1,3) + (2,4) + (3,4) 


(l,2,3) + 0,„,(l,3) + 0,„,(2,3) 
= (1,2, 3) + (2, 3, 4) 


(1,2,3) 









3.2. "Algebraic Thinning" 

Having obtained the thinned complex Mt„pK, we next construct a chain 
contraction {faig, Qaig, (paig) of the chain complex C{Mt„pK) to its homology. 
This step can be considered as a thinning at algebraic level (for this reason 
we call it "algebraic thinning"). We compute {faig, daig, (paig) interpreting the 
"incremental algorithm" j5| for computing homology groups in R^, in terms 
of chain contractions. As we will see later, the design of an algorithm for 
computing the cohomology ring of K, will be possible thanks to the infor- 
mation saved in the chain contraction of C{K) to its homology constructed 
before. 
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Let ((Ti, . . . , am) be a sorted set of all the simplices of K with the property 
that any subset {cii, . . . , cij}, i < m, is a sub complex of it. Algorithm 13.41 
computes a chain complex H with a set of generators h, and a chain con- 
traction {faig, daig, 4>aig) of (^(i^') to 71. Initially, k Is empty. In the ith step 
of the algorithm, the simplex cTj is added to the subcomplex {cri, . . . , cxi-i} 
and then, a homology class is created or destroyed. If faigd^Ci) = then 
(jj "creates" a homology class. Otherwise, (Tj "destroys" one homology class 
"involved" in the expression of faigd{ai). At the end of the algorithm, "H is 
a chain complex isomorphic to the homology of K. 
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Algorithm 3.4. Algebraic Thinning Algorithm 
Input: The sorted set {ai, . . . ,am)- 

Initially, faig{cr) := 0, 0„,g((7) := for each a e K; and /i := { }. 
For i = 1 to i = m do 
If falA(^^) = then 
h := hU{ai} 

Else take any one aj of faigd{ai) , then 

h:=h- {aj}, 

For k = 1 to A; = m do 

If aj appears in the expression of faig{o'k) then 

falgio'k) '■= falg{0-k) + falgd{ai) 
(t>alg{crk) (t>alg{cFk) + C^i + (t>algd{ai) 

End if. 
End for. 
End if. 
End for. 

For each cr e /i do 

9aig{cr) := (7 + 0„,g9((7). 
End for . 

Output: the chain contraction {faig, Qaig, <i>aig) of C{K) to H. 

The output of the algorithm allows us to determine both a representative 
cycle for each homology class and the homology class for each cycle. More- 
over, for any g-boundary a on K we can obtain a + l)-chain a' — (j)aig{ci) 
on K such that a — d{a'). 

Concerning to the complexity, suppose K has m simplices. In the ith 
step of the algorithm (1 < i < m), we have to evaluate dai. The number 
of simplices involved in dai is fewer or the same than the dimension of CTj 
which is at most 3. On the other hand, the number of elements involved 
in the formulae for faigdai and (j)aigdai is 0(3m) = 0{m). If dai 7^ 0, we 
have to update faigdak and (paigd^k for 1 < /c < m, so the total cost of these 
operations is 0{m^). Therefore, the total algorithm runs in time at most 

Let us observe that composing the chain contraction {ftop, Qtop, 4>top) of 
C{K) to C{Mt„pK), described in the previous subsection, with {faig, gaig, <l>aig) 
of C{MtopK) to H (isomorphic to H{K)), wc get a new chain contraction 

alg ftop") QtopQalg') ^top + gtop4>algftop) of C{K) to H. 
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Example 3.5. Let L be the simplicial complex showed in Figure O The 
intermediate stages of the algorithm are: 
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Finally, h = {(4), (1,2)}, g^,^{{4)) := (4) and g^,^{{l,2)) := (1,2) + 
(1,4) + (2, 3) + (2, 4). 

Summing up, the output is a chain contraction {faig, daig, (paig) of C{L) 
to the chain complex T-L^ (isomorphic to H{L)) with set of generators h = 
{(4), (1,2)}. In particular, we obtain that Ho{L) ~ Z/Z2, Hi{L) ~ Z/Z2 
and H2{L) = 0. 



3.3. Computing the Cohomology Ring 

After applying in order topological and algebraic thinning to the sim- 
plicial complex K, we are able to compute the multiplication table on the 
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Figure 6: The complex L (on the left) and the generators of %l (on the right). 

cohomology. Since the ground ring is a field, the homology and cohomology 
groups of K are always isomorphic. 

Let {f,g,4') be a chain contraction of C{K) to "H (where "H is a chain 
complex isomorphic to H{K) and to H*{K)), and let h = {ai, . . . ,ap} be 
a set of generators of "H obtained using the algorithms explained before. 
Then a*f : Cq{K) Z/Z2 (where a*{aj) = 1 H j = i and otherwise) 
is a representative cocycle of the cohomology class corresponding to a^, for 
I < i < p. 

Let a and (3 be two elements of h, then the cup product of the cohomology 
classes corresponding to a and (3 can be computed as follows: 

Algorithm 3.6. The cup product of two classes of cohomology. 

Input: the elements a and (3 

and the set of generators h = {ai, . . . ,ap} of Ti. 
Initially, := for 1 < k < p and cup := 0. 
For k = 1 to k = p do 

A,:=(a7-/37k(«fc). 
End for. 

cup := Yfk=l ^kOik- 

Output: the chain cup. 

Observe that the complexity of this algorithm for computing a ^ /3 is O(m^). 
Moreover, if we are interested in computing the cohomology ring of K., we 
have to apply the algorithm above to all the pairs {ai.aj), I < i < j < p 
(since the cup product is commutative, ai ^ aj = aj ai). Then, the 
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algorithm for computing the cohomology ring of K wiU run in time at most 
O(m^) if K has m simphces. 

Let us note that the cohomology ring of K is not suitable in general 
for topological classification tasks. This is due to the fact that determining 
whether two rings are isomorphic or not by means of their respective mul- 
tiplication tables is an extremely difficult computational question. In order 
to avoid this problem, we will put the information of the cup product table 
into a different form. 

If we restrict our interest in simplicial complexes are embedded in R^, 
observe that the possible non-trivial cup products are the ones j3 where 
both a and ^ are elements of h corresponding to cohomology classes of di- 
mension 1. 

In order to design a new algorithm for computing the cup product in a 
way that we can determine whether two cohomology rings are isomorphic or 
not by means of their respective multiplication tables, we need to define a new 
concept. Given a chain contraction (/, of C{K) to a chain complex l-L 
(isomorphic to H*{K)) with set of generators h, and a simplex a = {vq, ^1,^2) 
of dimension 2, suppose that {ai, . . . , ap} is the set of elements of h of di- 
mension 1, f{{vo,vi)) = J2iei (^i s-iid /((fi,f2)) = J2j£jC(j where / and J are 
subsets of the set {1, 2, . . . Define (/ /)(cr) = Y.i£i Y,j^j{oii-, oij). This 
definition can be extended to 2-chains by linearity. 

Algorithm 3.7. Cup Product Algorithm 

Input: A simplicial complex K 

and a chain contraction {f,g,4>) of C{K) 

to a chain complex (isomorphic to H*{K)) 

with set of generators h. 

Initially, q := the number of elements of h of dimension 2, 
6i := f or l<i<q and M ( ) . 

For i = 1 to i = q do 

hi := (/0/)5(«O- 
End for. 
M :=(&!,..., 6,). 
Output: The sorted set M. 

Let {ai, . . . , ap} be the set of elements of h of dimension 1 and . . . , Pq\ 
the ones of dimension 2. Each hi, 1 < ?' < is of the form -^}fc('^i) Qifc) 
where the sum is taken over the set {(j, k) : 1 < j,k < p) and A* = {a*f ^ 
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alf)g{f3i) could be or 1. Since the cup product is commutative, we have 
that A*;,. = XI.J. Therefore, the output of this algorithm can be put into a 
matrix form Ai of (cohomology classes 1 < i < q) x (pairs of cohomology 
clases {aj, ak), I < j < k < p). The column of Ai corresponding to the pair 
(ctj, ak), i < j < k < p, gives the value of the cup product aj ^ ak- This 
algorithm for computing the matrix A4 runs in time at most O(m^) if K has 
m simplices. 

From the diagonalization D of the matrix Ai, a first cohomology invari- 
ant HBi (K) for distinguishing non-homeomorphic simplicial complexes with 
isomorphic (co)homology groups appears. We define this cohomology num- 
ber in order to have a handy numerical tool for distinguishing 3D digital 
pictures. 

Definition 3.2. Given a simplicial complex K, the integer HBi(K) is a 
cohomology invariant defined as the rank of the matrix M. . 

4. A First Approach to the Digital Cohomology Ring 

Since an isomorphism of pictures on the 3D body centered cubic grid 
is equivalent to a simplicial homeomorphism of the corresponding simplicial 
representations, we are able to define the digital cohomology ring of / with 
coefficients in Z/Z2 as the cohomology ring of K[I) with coefficients in 7^/7/2. 
Moreover, the following definitions hold: 

Definition 4.1. Given a digital picture I — (V, 14, 14,S), the digital co- 
homology ring of I with coefficients in Z/Z2 is defined as the cohomology 
ring of K{I) with coefficients in Z/Z2. The cohomology invariant HBi{I) 
is defined as HBi{K{I)). 

In the previous sections, we have showed that it is possible to com- 
pute the digital cohomology ring of / with coefficients in Z/Z2. The steps 
of the method are: first, we construct the simphcial complex K{I). Sec- 
ond, we topologically thin K{I), obtaining a smaller simplicial complex 
M,,pK{I) and a chain contraction (/^^p, g^.^, (p^op) of C{K{I)) to C{Mt,^K{I)). 
Third, we compute "H which is isomorphic to H{I) and a chain contraction 
(/oig, Qaig, 4>aig) of C {M^^pK {!)) to H. FouT , wc calcuktc the cohomology ring 
of / via the cohomology ring of Mt„pK{I) and the invariant HBi{I) via 
H Bi{Mt„pK (I)) , using the chain contractions constructed before. All the 
information obtained in this way is useful for topologically classifying and 
distinguishing binary 3D digital pictures. 
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4.1. Some Examples 

In order to show examples of the computation and visuahzation of the 
cohomology ring of simple pictures, we expose a small prototype called Ed- 
itCup. We use a free program for building 3D worlds. In our case, a world 
is a particular 3D simplicial complex K representing a digital picture I con- 
sidering the 14-adjacency. A way for distinguishing the different maximal 
simplices of a simplicial representation is by using different colours: red for 
tetrahedra, green for triangles, blue for edges, and black for vertices. 

All the computations are done considering Z/Z2 as the ground ring. For 
visualizing (co) chains, the simplices on which a given (co) chain is non-null, 
are lighted in a different color. On the other hand, the "visualization" of any 
(co)homology class on K is given by lighting the simplices of K on which the 
representative cochain of this class is non-null. Moreover, the "visualization" 
of any (co) homology class on the original 3D digital binary- valued picture / 
could be given by lighting the points of / such that the corresponding vertices 
span simplices on which the representative cochain of this class is non-null. 

The DPS used in these examples, that we call (14, 14)-DPS, is (Z^ 14, 14), 
in which the underlying grid is the set of points with integer coordinates 
in the Euclidean 3-space and the 14-neighbours of a grid point (black 
or white) with integer coordinates {x,y,z) are: (x ± l,y,z), {x,y ± 1,^), 
{x,y,z±l), {x + l,y-l,z), (x - 1, ?/ + 1, z), {x + l,y,z-l), (x - 1, z + 1), 
{x,y + l,z-l), {x,y-l,z + l), {x + l,y + 1, z - 1), {x - l,y - 1, z + 1) (see 
Figure [7]). The (14, 14)-DPS and the BCC grid are isomorphic DPSs: each 
grid point (x, y, z) of the (14, 14)-DPS can be associated to a point (a, 6, c) 
via the formula: (a, 6, c) = (x + y — 2z, —x + y, — x — y)). 

Let us consider now the following pictures: a torus I = (Z^, 14, 14, 5^^) 
and a wedge of two topological circles and a topological 2-sphere J = 
(Z^, 14, 14, Bj) (see Figure[8]). In the volumetric representation of the picture 
/ (resp. J), we use voxels with centres the points Bj (resp. Bj). It is clear 
that the (co) homology groups of I are isomorphic to those of J. They are 
Z/Z2, Z/Z2 0Z/Z2 and Z/Z2 of dimension 0, 1 and 2, respectively. So, the 
(co)homology information is not enough for topological distinguishing both 
pictures. 

Let K{I) and K{J) be the simplicial representations of / and J respec- 
tively (see Figure In order to compare the cohomology ring of both 
pictures, the first step is the computation of chain contractions {fi,gi,4>i) 
of C{K{I)) to "Hi and {fj,gj,4>j) of C{K{J)) to Tij using the topological 
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Figure 7: The 14-neighbours of a grid point p of the (14,14)-DPS. 




Figure 8: The volumetric representation of the pictures / (on the left) and J (on the 
right). 

and algebraic thinning algorithms explained before, where T-Lx (resp. T-Lj) is 
a chain complex isomorphic to the (co)homology of / (resp. J). 

Let us denote by ai and 02 (resp. a'^ and aj) the generators of l-ix 
(resp. l-ij) of dimension 1 and 03 (resp. 0^3) the generator of l-Lx (resp. 
T-Lj) of dimension 2. Let us also denote by the representative cycles of 
the generators of H{I) (that is, = gj{ai))] and by the same of H{J). 
We visualize these cycles on K{I) and K{J) in Figure [TOl In Figure [TT| the 
representative cocycles bi (resp. 6-) obtained via the formula bi = a*fi (resp. 
b'i = Ci'* fj) of the generators of H*{I) (resp. H*[J)) are shown. Recall that 
we do it by lighting the simplices on which the cochains is non-null. 

The output of Algorithm 13 . 71 for "H/ and are Mj = ((ai, 02) + (^2, cn)) 
and Mj = (0), respectively. The matrices corresponding to the cohomology 
rings of the pictures / and J are: 
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Figure 9: The simplicial complexes K{I) (on the left) and K{J) (on the right). 




Figure 10: The cycles oi, 02 and a'^, al^ (in yellow); and 03 and (in green). 



Figure 11: The cocycles 61, hi and h'^, h'2 (in yellow); and 63 and 63 (in green). 
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(ai,ai) ioLx.OL'i) (a2,a2) 
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(a'i,a'i) 


(a;,a^) 


(a2,«2) 


as 


1 


"3 












Therefore, EBx{l) = 1 and HBi{J) = 0. We conclude that K{I) and 
K[J) are not homeomorphic (more precisely, we conclude that they are not 
homotopy equivalent), then I and J are not isomorphic. 

Let us expose another example (see Figure [12]): the picture A is a wedge 
of two torus; the picture B consists in a wedge of a sphere and a genus-2 torus 
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Figure 12: The pictures A (on the left) and B (on the right). 




Figure 13: The simphcial complexes K{A) (on the left) and K{B) (on the right). 

(a sphere with two handles and two holes). Both pictures have 1 connected 
component, 4 holes and 2 cavities. 

The simplicial representations of A and B, K{A) and K{B) are showed 
in Figure [T3l In Figure [H] (resp. Figure [T5|) . the representative cycles and 
cocycles of the generators of the (co)homology of A (resp. B) are showed. 

Let us denote by (resp. a'^), i = 1,2,3,4, the generators of Ha (resp. 
Tig) of dimension 1; and by f3i (resp. (3^), i = 1,2, the generators of Ha (resp. 
Hg) of dimension 2. All of them are obtained using the algebraic thinning 
algorithm explained in the previous section. The output of Algorithm 13 . 71 for 
Ha and Hg is Ma = {{ai, a2) + {ai, a3) + {a2, ai) + (a3, ai), (as, a4) + (a4, 03)) 
and Mg = (0, (ai a^) + (0:2, as) + (as, 0^2) + («45 Q^i)) • Therefore, the matrices 
corresponding to the cohomology rings of A and B are: 
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Figure 14: On the left (resp. on the right), the representative cycles (resp. cocycles) of 
the generators of H{A) (resp. of H*{A)). 



p5¥ 

Figure 15: On the left (resp. on the right), the representative cycles (resp. cocycles) of 
the generators of H{B) (resp. of H*{B)). 
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where (z, j) represents the pair {ai^aj) (resp. {a'^^a'^)). 

We conclude that HBi{A) = 2 and HBi{B) = 1, and then K{A) and 
K{B) are not homeomorphic (more precisely, we conclude that they are not 
homotopy equivalent), therefore A and B are not isomorphic. 



5. Conclusions and Future Work 

We have seen that there is a true algorithm for computing the digital co- 
homology ring (with coefficients in Z/Z2) of a 3D binary picture on the BCC 
grid. It is also possible to compute the digital cohomology ring of / with coef- 
ficients in any commutative ring G, thanks to the fact the simplicial complex 




22 



K{I) is embedded in and, consequently, it have torsion-free homology. 
We deal here with Z/Z2 coefficients, in order to simplify and avoiding signs 
in the explanation of our algorithmic formulation, to give an easy geomet- 
ric interpretation of "digital" cohomology classes and to work with binary 
arithmetic. Moreover, there is no problem to define the cohomology ring of 
/ with coefficients in a commutative ring G as the cohomology ring of K{I) 
with coefficients in G; and the cohomology invariant HBi{I\ G) with coeffi- 
cients in G as HBi{K{I);G). On the other hand, since HBi{K{I); G) can 
be obtained from the ffist homology group of the reduced bar construction 
B{C*iK{I);G)) Q associated to the cochain complex C*{K{iy,G) with 
coefficients in G, we will confine ourselves to say that the rest of homology 
groups of this last algebraic object give rise to more complicated cohomology 
invariants for a digital binary-valued picture. 

In this paper, we talk about topological and algebraic thinning. Concern- 
ing the ffist one, we do not use here well-known direct (without passing to 
simplicial framework) topological thinnings of digital binary-valued pictures 
because we are interested in constructing chain contractions which allow us 
to obtain cohomology results. Concerning the second one, the idea of com- 
puting a chain contraction of a chain complex to its homology has also been 
used in [3] for computing primary and secondary cohomology operations. 

Another important question is to try to improve the complexity of the 
algorithm computing digital cohomology ring on the BCC grid detailed in 
this paper. We do not take advantage here neither of the particular sim- 
plicial structure of the simplicial complex K{I) (determined by the BCC 
grid) associated to /, nor of representing in a compressed form (without loss 
of information) the 3D digital picture (for example, in an octree format). 
To obtain positive results in these directions and to eliminate from our al- 
gorithm the intermediary simplicial objects will allow us to specify a more 
refined algorithm computing digital cohomology on the BCC grid. 

Finally, another important question that it is necessary to deal with in a 
near future is to try to generalize this work to other natural digital picture 
spaces. 

References 

[1] Ayala R., Dommguez E., Frances A.R., Quintero A.: Homotopy in Dig- 
ital Spaces. DGCl 2000. LNCS, Springer-Verlag 1953 (2000) 3-14 

[2] AlexandrofF P., Hopf H.: Topologie I. Springer, Berlin 1935 



23 



[3] Berrio J.M., Gonzalez-Diaz R., Leal F., Lopez M.M., Real P.: Visual- 
izing Cohomology Aspects of 3D Objects. Proc. of the 6th Asian Tech. 
Conf. in Math. (2001) 459-468 

[4] Bjorner A.: Topological Methods. Handbook on Combinatorics. Elsevier 
Sci. 2 (1995) 1819-1872 



24 



[5] Delfinado C.J. A., Edclsbrunner H.: An Incremental Algorithm for Betti 
Numbers of Simplicial Complexes on the 3-Sphere. Comput. Aided 
Geom. Design 12 (1995) 771-784 

[6] Forman R.: Combinatorial Differential Topology and Geometry. New 
Perspective in Geom. Combinatorics. MSRI Pubhc. 8 (1999) 177-206 

[7] Gonzalez-Diaz R., Real R: Computation of Cohomology Operations on 
Finite Simplicial Complexes. Homology, Homotopy and Applications 5 
(2) (2003) 83-93 

[8] Gonzalez-Diaz R., Real P.: Towards Digital Cohomology. DGCI 2003, 
LNCS, Springer 2886 (2003) 92-101 

[9] Kenmochi Y., Imiya A.: Polyhedral Set Operations for 3D Discrete 
Object Deformation. Scale-Space 2001, LNCS, Springer-Verlag and 
lEEE/CS 2106 (2001) 417-425 

[10] Khahmsky E.D., Kopperman R.D., Meyer P.R.: Computer Graphics 
and Connected Topologies on Finite Ordered Sets. Topology and Appl. 
36 (1990) 1-17 

[11] Kong T.Y.: A Digital Fundamental Group. Computer Graphics 13 
(1989) 159-166 

[12] Kong T.Y., Roscoe A.W., Rosenfeld A.: Concepts of Digital Topology. 
Topology and its Applications 46 (1992) 219-262 

[13] Kovalevsky V.A.: Discrete Topology and Contour Definition. Pattern 
Recognition Letter 2 (1984) 281-288 

[14] MacLane S.: Homology. Classic in Math., Springer-Verlag 1995 

[15] Munkres J.R.: Elements of Algebraic Topology. Addison- Wesley Co. 
1984 

[16] Rosenfeld A.: 3D Digital Topology. Inform, and Control 50 (1981) 119- 
127 

[17] Theussl T., MoUer T., GroUer M.E.: Optimal Regular Volume Samphng. 
Proc. IEEE Conference on Visualization (2001). 



25 



